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Abstract 



We discuss to which extent the present experiments of direct search for 
WIMPs, when interpreted in terms of rehc neutrahnos, probe interesting 
regions of the supersymmetric parameter space, which are also being pro- 
gressively explored at accelerators. Our analysis is performed in a number of 
different supersymmetric schemes. We derive the relevant neutralino cosmo- 
logical properties, locally and on the average in the universe. We prove that 
part of the susy configurations probed by current WIMP experiments entail 
relic neutrahnos of cosmological interest. The main astrophysical and particle 
physics uncertainties, relevant for a proper comparison between theory and 
experimental data, are stressed and taken into account. 
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I. INTRODUCTION 



As was first noticed in Ref . , in the last few years tlie experiments of direct searcfi for 
Weakly Interacting Massive Particles (WIMP) p[ have already reached a sensitivity which 
allows the exploration of regions of the supersymmetric parameter space, which are also 
progressively investigated at accelerators. This property is manifest, when the experimental 
results are interpreted in terms of relic neutralinos [||]. 

The probing of the susy parameter space by WIMP direct searches is even more sizeable 
at present, with the sensitivities of experiments PJ^; a comparative discussion of the experi- 
mental features and implications of the DAMA ^ and CDMS ^ experiments may be found 
in Ref. 0. Detailed studies of the possible interpretation of the annual-modulation effect 
in terms of relic neutralinos have been reported in 1^-0. Comparisons of the experimental 
data of Ref. with susy calculations have also been performed in Refs. [Il0|-|15|. 

In the present paper we intend to clarify the actual capability of WIMP direct searches by 
exploring in a sistematic way different realizations of minimal supersymmetric models and 
showing their intrinsic differences in the prediction of neutralino rates and relic abundance. 
Specifically, we will consider two different implementations of a supergravity scheme with 
parameters defined at the grand unification scale and an effective supersymmetric model 
defined at the electroweak scale. Our analyses are performed in the light of the follow- 
ing relevant points: i) current uncertainties in astrophysical properties, ii) uncertainties 
in hadronic quantities, iii) new bounds from LEP searches for Higgs and supersymmetric 
particles, iv) updated determinations of cosmological parameters. 

Let us start by recalling that the determination of the sensitivity range of an experiment 
of WIMP direct search in terms of the WIMP mass and of the WIMP-nucleon cross section 
rests on a number of crucial assumptions, since it depends both on the distribution function 
of the WIMPs in the halo and on the nature of the relic particle. 

A WIMP direct experiment provides a measurement (or an upper bound) of the differ- 
ential event rate 

§- = NrJ^fdvf(v)v-^(v.E^) (1) 

where is the number of the target nuclei per unit of mass, m\Y is the WIMP mass, pw is 
the local WIMP matter density, v and f{v) denote the WIMP velocity and velocity distri- 
bution function in the Earth frame {v = \v\) and da/dEn is the WIMP-nucleus differential 
cross section. The nuclear recoil energy is given by Er = rn^^^v^^l — cos9*)/mN, where 9* 
is the scattering angle in the WIMP-nucleus center-of-mass frame, mM is the nuclear mass 
and rrij-ed is the WIMP-nucleus reduced mass. Eq.(0) refers to the case of a monoatomic 
detector, like the Ge detectors. Its generalization to more general situations, like for instance 
the case of Nal, is straightforward. In what follows pw will be factorized in terms of the 
local value for the total non-baryonic dark matter density pi and of the fractional amount 
of density, ^, contributed by the candidate WIMP, i.e. pw = ^ ■ Pi- For pi we use the range 
0.2 GeV cm"'^ ^ Pi ^ 0.7 GeV cm~^, where the upper side of the range takes into account 
the possibility that the matter density distribution is not spherical, but is described by an 
oblate spheroidal distribution |T6|,0. 
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The WIMP-nucleus differential cross section may conveniently be split into a coherent 
part and a spin-dependent one 



da 




(2) 



SD 



whose generic features are discussed in the seminal paper of Ref . . To compare theoretical 



expectations with experimental data, and experimental data of different detectors among 
themselves, it is useful to convert the WIMP-nucleus cross-section into a WIMP-nucleon 
cross section. This procedure is feasible independently of the nuclear model and of the 
specific nature of the WIMP only under the hypothesis that the coherent cross-section is 
dominant and the WIMP couples equally to protons and neutrons (at least approximately) 
||l||. Under this assumption, the WIMP-nucleus cross section may be expressed in terms of 
a WIMP-nucleon scalar cross section cr^caiar'"'^'' 



da ^ ( da \ ^ F'^{q) ( I + mw/mp^ 



4 2 (nucleoli) /„n 
^ "^scalar ) K'-') 



where irip and mi^ are the proton and nucleus mass, A is the nuclear mass number, E^°'^ 
is the maximal recoil energy and F{q) is the nuclear form factor for coherent interactions. 



This form factor is usually parametrized in the Helm form [|1^]; however, precise evaluations 
of the event rates may require specific nuclear calculations for each target nucleus. In the 
rest of this paper we assume that the WIMP interaction with the nuclei of the detector is 
dominated by coherent effects, so that a WIMP-nucleon scalar cross section may be derived 
from the WIMP-nucleus cross section by use of Eq. (^. 

Now, coming back to the general expression in Eq. (|ip, we stress that extracting an 
information about the WIMP-nucleus cross sections from the experimental data requires the 
use of a specific expression for the velocity distribution function f[v) (notice that in writing 
Eq. (0) we have already made the assumption that the WIMP phase-space distribution 
function may be factorized as p ■ f{v), and this is certainly not the most general case [|l6|] ). 
The usual choice for f{v) is the isotropic Maxwell-Boltzmann distribution in the galactic 
rest frame, as derived from the isothermal-sphere model. 

However, recent investigations have shown that deviations from this standard scheme, 



either due to a bulk rotation of the dark halo [^,0 or to an asymmetry in the WIMP ve 



locity distribution P^|-p^, infiuence the determination of the WIMP-nucleus cross sections 
from the experimental data in a sizeable way. In Ref. also triaxial matter distributions 



are considered; in the present paper deviation from sphericity in the WIMP matter distribu- 
tions are taken into account only through the physical range allowed for the value of pi (see 
our previous comment on pw after Eq. (|l|)). In a typical plot, where the WIMP-nucleus 
cross section is given in terms of the WIMP mass, the effect introduced by the mentioned 
deviations from the Maxwell-Boltzmann distribution is generically to elongate the contours 
towards larger values of mw This is for instance the case for the the annual-modulation 
region of the DAMA Collaboration 0. In Fig. 3 of Ref. it is shown that, by imple- 
menting the dark halo with a bulk rotation according to the treatment in Ref. |]2T|, the 
annual-modulation region moves towards larger values of the WIMP mass, with an elonga- 
tion which brings the right-hand extreme from the value of ~ 150 GeV to ~ 200 GeV. A 
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similar effect is obtained by introducing an asymmetry in the WIMP velocity distribution 



f{v)- Fig. 4 of Ref. |24| illustrates this point. Notice that this asymmetry effect also pushes 
somewhat downwards the annual-modulation region. We emphasize that all these effects 
are extremely important, when experimental results of WIMP direct detection are being 
compared with theoretical models for specific candidates. This point has been overlooked 
in most analyses in terms of relic neutralinos p5 |. 



In the present paper we focus our analysis to the WIMP mass range which, in the light of 
the present experimental data and of the previous considerations on the astrophysical 
uncertainties, appears particularly appealing: 

40 GeV <mw< 200 GeV. (4) 

Let us notice that the mass range of Eq. is quite appropriate for neutralinos. Actually, 
the lower extreme is indicative of the LEP lower bound on the neutralino mass (in the 
calculations performed in the present work the actual lower bound for m^, dependent on 
the other susy parameters, is employed, according to the constraints given in [Q). As for 
the upper extreme, we notice that, though a generic range for might extend up to about 
1 TeV, requirements of no excessive fine-tuning would actually favour an upper bound 
of order 200 GeV, in accordance with Eq. (^. 

In what follows we will discuss the discovery potential of WIMP direct searches for 
WIMPs in the mass range of Eq.(^. Particular attention will be paid to capabilities of the 
present experiments; their sensitivity range, in case of WIMPs whose coherent interactions 
with ordinary matter are dominant over the the spin-dependent ones, may be stated, in 
terms of the quantity ^o■^c^lar°''^ as 

4 ■ 10-1° nbarn < e^S^"^ < 2 • 10"' nbarn. (5) 

We will hereafter refer to region R as the one in the space mw —^o'^^^ilr"^^ which is defined 



by Eqs. (§-§[). The region R represents the sensitivity region already under exploration with 
present detectors. 

Our analysis, based on an interpretation of experimental data in terms of relic neutrali- 
nos, will show by how much the WIMP direct searches probe the supersymmetric parameter 
space. We remark that, in the case of neutralinos, the assumption about the dominance 
of the coherent cross section over the spin-dependent one is, in general, largely satisfied, 
except for values of cr^^j^^^^""^ which are far below the present experimental reach [|T|. 

The present analysis will be performed in the framework of various schemes, from those 
based on universal or non-universal supergravity, with susy parameters defined at the grand 
unification scale, to an effective supersymmetric model defined at the Electro- Weak (EW) 
scale. This is discussed in Sect. II, where we also specify the values employed here for the 
Higgs-quark-quark and the neutralino-quark-squark couplings. These quantities are sub- 
ject to sizeable uncertainties, as was stressed in Ref. p8|, which triggered a reconsideration 
of this important point in a number of subsequent papers |T^,|T^,^. 

The most important properties to be established for the relic neutralinos, which are en- 
tailed in the exploration by WIMP direct searches, concern their cosmological properties. 
Here, we perform a general analysis which is not limited to a restricted range of the cos- 
mological matter abundance, flmh"^ i^m is the matter cosmological density divided by the 
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critical density and h is the Hubble constant in units of 100 km s~^ Mpc~^) . Instead, we 
derive the average and local cosmological properties of the susy configurations from experi- 
mental determinations of c^gcalar'°"^ without any a priori requested range on Vtmh"^- On the 
basis of the results of our evaluations in the various supersymmetry models, we discuss when 
the relic neutralino does or does not saturate the expected amount of the local and of the 
average amount of total dark matter. Our results and conclusions are presented in Sects. 
Ill and IV, respectively. 



II. SUPERSYMMETRIC MODELS 

The calculations presented in this paper are based on the Minimal Supersymmetric 
extension of the Standard Model (MSSM), in a variety of different schemes. The essential 
elements of the MSSM are described by a Yang-Mills Lagrangian, the superpotential, which 
contains all the Yukawa interactions between the standard and supersymmetric fields, and 
by the soft-breaking Lagrangian, which models the breaking of supersymmetry. To fix the 
notations, we write down explicitly the soft supersymmetry breaking terms 



2| 



+ { [A[,h[bLaH,Rb + AiXbQaHibb + AXbQaH2Ub + h.c] - BfiHiH2 + h.c.} 
+ 5]M,(AiA, + AiA,) (6) 

i 

where the 0j are the scalar fields, the Aj are the gaugino fields. Hi and H2 are the two Higgs 
fields, Q and L are the doublet squark and slepton fields, respectively, and U, D and R 
denote the 5'f/(2)-singlet fields for the up-squarks, down-squarks and sleptons. In Eq.(|), 
nti and Mi are the mass parameters of the scalar and gaugino fields, respectively, and A 
and B denote trilinear and bilinear supersymmetry breaking parameters, respectively. The 
Yukawa interactions are described by the parameters which are related to the masses of the 
standard fermions by the usual expressions, e.g., rrit = h^V2, rnb = h'^Vi, where Vi =< Hi >. 

Implementation of this model within a supergravity scheme leads naturally to a set of 
unification assumptions at a Grand Unification (GUT) scale, Mqut'- 

i) Unification of the gaugino masses: Mi^McuT) = ^1/2, 

ii) Universality of the scalar masses with a common mass denoted by mg: mi{MGUT) 

= mo, 

iii) Universality of the trilinear scalar couplings: A^Mqut) = A'^{Mgut) = A'^{Mgut) 

= Aonio. 

This scheme will be denoted here as universal SUGRA (or simply SUGRA). The relevant 
parameters of the model at the electro-weak (EW) scale are obtained from their correspond- 
ing values at the Mqut scale by running these down according to the renormalization group 
equations (RGE). By requiring that the electroweak symmetry breaking is induced radia- 
tively by the soft supersymmetry breaking, one finally reduces the model parameters to five: 
mii2, ruQ, Aq, tan I3{= V2/V1) and sign fi. In the present paper, these parameters are varied 
in the following ranges: 50 GeV < mi/2 < 1 TeV, mo < 1 TeV, —3<Ao< +3, 1 < 
tan /3 < 50. Notice that a common upper extreme for the mass parameters has been used, 
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and generically set at the value of 1 TeV, as a typical scale beyond which the main attractive 
features of supersymmetry fade away. However, fine-tuning arguments actually set different 
bounds for tjiq and mi/2 (in universal SUGRA and in nuSUGRA) [^: mi/2 < hundreds 
of GeV, whereas mo < 2-3 TeV. In the present paper we did not look specifically into the 
mo ~ 2 — 3 TeV window; in Ref. |^ phenomenology of relic neutralinos in this large mo 
regime has been analyzed |^ 



Models with unification conditions at the GUT scale represent an appealing scenario; 
however, some of the assumptions listed above, particularly ii) and iii), are not very solid, 
since, as was already emphasized some time ago universality might occur at a scale 
higher than Mqut ~ 10^^ GeV, e.g., at the Planck scale. More recently, the possibility that 
the initial scale for the RGE running, M/, might be smaller than Mqut ~ 10^^ has been 
raised |T^,|33|, on the basis of a number of string models (see for instance the references 



quoted in [|T4|)- In Ref. ||T^ it is stressed that Mj might be anywhere between the EW scale 
and the Planck scale, with significant consequences for the size of the neutralino-nucleon 
cross section. 

An empirical way of taking into account the uncertainty in Mj consists in allowing 
deviations in the unification conditions at Mqut- For instance, deviations from universality 
in the scalar masses at Mquti which split Mh^ from M^.^ may be parametrized as 

Ml^{MGUT) = ml{l + 5,) (7) 

This is the case of non-universal SUGRA (nuSUGRA) that we considered in Refs. [|7|,|27|, 
and that we analyse again in this paper. Here the parameters 5j which quantify the depar- 
ture from universality for the M^. will be varied in the range (-2, +2). Deviations from 
universality in the Higgs masses have recently been considered also in Ref. |T5[. Further 



extensions of deviations from universality in SUGRA models which include squark and/or 
gaugino masses are discussed, for instance, in IT^JSS 



The large uncertainties involved in the choice of the scale M/ make the use of SUGRA 
schemes rather problematic and unpractical: the originally appealing feature of a univer- 
sal SUGRA with few parameters fails, because of the need to take into consideration the 
variability of Mj or, alternatively, to add new parameters which quantify the various devi- 
ation effects from universality at the GUT scale. It appears more convenient to work with 
a phenomenological susy model whose parameters are defined directly at the electroweak 
scale. We denote here this effective scheme of MSSM by effMSSM. This provides, at the EW 
scale, a model, defined in terms of a minimum number of parameters: only those necessary 
to shape the essentials of the theoretical structure of an MSSM, and of its particle content. 
Once all experimental and theoretical constraints are implemented in this effMSSM model, 
one may investigate its compatibility with SUGRA schemes at the desired M/. 

In the effMSSM scheme we consider here, we impose a set of assumptions at the elec- 
troweak scale: a) all trilinear parameters are set to zero except those of the third family, 
which are unified to a common value A\ b) all squark soft-mass parameters are taken de- 
generate: rrtq. = mq] c) all slepton soft-mass parameters are taken degenerate: m^- = mj; 
d) the U{1) and SU{2) gaugino masses. Mi and M2, are assumed to be linked by the usual 
relation Mi = (5/3) tan^ ^^^^2 (this is the only GUT-induced relation we are using, since 
gaugino mass unification appears to be better motivated than scalar masses universality). 
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As a consequence, the supersymmetric parameter space consists of seven independent pa- 
rameters. We choose them to be: M2, tan/?, m^, m^, m^-, A and vary these parameters in 
the following ranges: 50 GeV < M2 < 1 TeV, 50 GeV < < 1 TeV, 80 GeV < < 
1 TeV, 100 GeV < m^, m^- < 1 TeV, -3 < A < +3, 1 < tan /? < 50 (m^ is the mass of the 
CP-odd neutral Higgs boson). 

The effMSSM scheme proves very manageable for the susy phenomenology at the EW 
scale; as such, it has been frequently used in the literature in connection with relic neu- 
tralinos (often with the further assumption of slepton-squark mass degeneracy: rriq = mj) 
PJ9| jTI|J5^ -|5B| . Notice that we are not assuming here slepton-squark mass degeneracy. In 
the scatter plots given in this paper only configurations with rriq > mj are shown. This mass 
hierarchy is reminiscent of what is usually obtained in SUGRA schemes, although in our 
effMSSM it is not necessarily so. It is worth reporting that some configurations with inverse 
hierarchy rriq < m^- produce some increase in ^cr^caiar'"'^'' ^^"^ '^x values (see the discussion 
after Figs. 5 in Sect. III). 

We recall that even much larger extensions of the supersymmetric models could be en- 



visaged: for instance, non-unification of the gaugino masses |29,3^, and schemes with CP- 



violating phases ||38|. Here we limit our considerations to the schemes previously defined: 
universal SUGRA, nuSUGRA, effMSSM. 

The neutralino is defined as the lowest-mass linear superposition of photino (7), zino 
(Z) and the two higgsino states {H^, H2): x = '^i7 + '^2^ + a-iHl + aa^H^. Hereafter, the 
nature of the neutralino is classified in terms of a parameter P, defined a.s P = a\ + a\. The 
neutralino is called a gaugino when P > 0.9, a higgsino when P < 0.1, mixed otherwise. 

For more details concerning theoretical aspects involved in our calculations and the way 
in which the experimental constraints due to 6 — s + 7 is implemented we refer to Refs. 
I^Q]. Accelerators data on supersymmetric and Higgs boson searches (CERN e'^e~ collider 
LEP2 and Collider Detector CDF at Fermilab) provide now rather stringent bounds on 
supersymmetric parameters. CDF bounds are taken from [R^. The new LEP2 bounds are 



taken from p6| , |40[| ; these constrain the configurations of relevance for relic neutralinos more 



severely as compared, for instance, with those considered in Ref. [^. 

The results for the neutralino relic abundance have been obtained with the procedure 
indicated in Ref. [^. The neutralino-nucleon cross section has been calculated with the 



formulae reported in Ref. PJ28[]. As discussed in the introduction, this cross section suffers 
from significant uncertainties in the size of Higgs-quark-quark and squark-quark-neutralino 
couplings. In this paper we use for these quantities what we have defined as set 1 and set 
2 in Ref. to which we refer for details. Here we only report the values of the quantities 



rriq < qq> for the two sets: 

Set 1: mi <ll> = 23 MeV, m, < ss > = 215 MeV, rrih < hh > = 50 MeV. (8) 



Set 2:mi<U>= 30 MeV, m, < ss > = 435 MeV, < hh > = 33 MeV. (9) 

In Eqs. (^ ^ / stands for light quarks, s is the strange quark and h = c,b,t denotes heavy 
quarks. For the light quarks, we have defined mi <ll > = \[mu < uu > +md < dd >]. Set 
1 and set 2 bracket, only partially, the present uncertainties. In Ref. we also considered 
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the consequences of using a more extreme set of values. It is worth recalhng that the quantity 
rris < ss > is crucial in establishing the size of CTgcaiar""'' [H]- 

The results shown in the next section are obtained with the same numerical codes em- 
ployed in our previous papers [0-P,|28l 



but take into account all new accelerator data. 



III. RESULTS 

We turn now to the presentation of our results. In Figs, la-c we give the scatter plots 
for cTgcaiar""^ vcrsus ^x^'^ ^^^^ three different schemes: universal SUGRA, non-universal 
SUGRA and effMSSM. For the SUGRA schemes we only display the results corresponding 
to positive values of /x, since, for negative values, the constraint on 6 — >• s + 7 implies a 
large suppression of Cg^^i^^^"'^'* . The two horizontal lines bracket the sensitivity region defined 
by Eq. (^. The two vertical lines denote a favorite range for Qmh'^, 0.05 < Vtmh'^ < 0.3, 
as derived from a host of observational data. According to the most recent determinations 
43| , the lower bound on VL^h? is approaching the value 0.08. However, due to the still 



unsettled situation as regards determinations of the matter density in the universe and of 
the Hubble constant, a conservative attitude seems advisable. Anyway, we stress that in 
the present paper we are not restricting ourselves to any particular interval of ^mh"^- Only 
some features of Figs. 5 depend on the actual value employed for the minimum amount of 
matter necessary to reproduce the halo properties correctly. 

Figs, la-c provide a first relevant result of our analysis: the present experimental sen- 
sitivity in WIMP direct searches allows the exploration of supersymmetric configurations 
of cosmological interest, also in the constrained SUGRA scheme. It is remarkable that the 
upper frontier of the scatter plots is not significantly different in the three different models, 
although the region of experimental sensitivity and cosmological interest is covered with an 
increasingly larger variety of supersymmetric configurations as one moves from SUGRA to 
nuSUGRA and to effMSSM. This latter fact is expected from the intrinsic features of the 
various schemes. This point will be further discussed later on, in connection with Figs. 4. 
Fig. 2 shows what is the effect of using set 2 instead of set 1 for the quantities rriq < gg >'s 
in effMSSM. 

Once a measurement of the quantity Px ' '^scalar'""'' is performed, values for the local density 



Px versus the relic abundance f^^/i^ may be deduced by proceeding in the following way p8| 

H \ . 1 J. 1 r (nucleoli) n / (nucleon) i r (nucleon) n i , ,1 

1) Px IS evaluated as [p^ ■ <aiar Uxpt / <aiar . where [p^ ■ a^'^aiar Uxpt denotes the 



experimental value, and crg™!^^^"'^'' is calculated as indicated above; 2) to each value of p^ 
one associates the corresponding calculated value of fi^/i^. The scatter plot in Fig. 3 is 
derived from the lowest value of the annual-modulation region of Ref. [^, [p^/(0.3 GeV 

cm~^) ■o'l^ailv°^^]expt = 1 " 10~^ ubam, and by taking in the range of Eq. (H). This plot, 
obtained in case of effMSSM, shows that the most interesting region, i. e. the one with 
0.2 GeV cm~'^ < Px — ^-"^ GeV cm~^ and 0.05 < Qmh^ < 0.3 (cross-hatched region in the 
figure), is covered by susy configurations probed by the WIMP direct detection. 

Let us examine the various sectors of Fig. 3. Configurations above the upper horizontal 
line are incompatible with the upper limit on the local density of dark matter in our Galaxy 
and must be disregarded. Configurations above the upper slanted dot-dashed line and below 
the upper horizontal solid line would imply a stronger clustering of neutralinos in our halo 
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as compared to their average distribution in the Universe. This situation may be considered 
unhkely, since in this case neutrahnos could fulfill the experimental range for p^, but they 
would contribute only a small fraction to the cosmological cold dark matter content. For 
configurations which fall inside the band delimited by the slanted dot-dashed lines and 
simply-hatched in the figure, the neutralino would provide only a fraction of the cold dark 
matter at the level of local density and of the average relic abundance, a situation which 
would be possible, for instance, if the neutralino is not the unique cold dark matter particle 
component. To neutralinos belonging to these configurations one should assign a resettled 
local density Px = Pi ^ ^x^'^ / (^mh'^)min, where {Qmh'^)min is the minimum value of flmh^ 
compatible with halo properties. 

It is interesting to analyze the properties pertaining to the supersymmetric configurations 
which stay inside the favored region 0.2 GeV cm~^ < Px — ^-"^ GeV cm~^ and 0.05 < flmh"^ < 
0.3 and to those which stay inside the corridor where rescaling applies {i.e. the corridor 
between the two slanted dot-dashed lines). In Figs. 4a-c we give the scatter plots for all 
these configurations in the plane rrih — tan f3 {rrih being the mass of the lightest CP-even 
neutral Higgs boson). 

Let us make a few comments about these results. First, we note that a feature that was 
already pointed out in Ref. ^ is recovered: in SUGRA (Fig. 4a) only high values of tan (3, 
tan/3 > 40, are involved in present direct detection experiments. Similar conclusions are 
also reached in more recent papers 0,0] . The occurrence of the lower bound tan /5 > 40 in 
the SUGRA scheme is a consequence of RGE evolutions of the parameters and the nature 
of radiative electroweak symmetry breaking, which induce strong correlations among the 
parameters at the low energy scale [^. In this class of models, couplings of the light Higgs 
boson h to the s quark can be enhanced only for large values of tan/?. As a consequence, 
the neutralino-nucleon cross section can be substantially large only when tan /3 > 40. Con- 
figurations displayed in Fig. 4a entail relatively light pseudoscalar Higgs A, niA < 200 GeV, 
'^o ^ 350 GeV, and lightest stop and sbottom masses larger than about 460 GeV. 

These features are somewhat (even though not completely) relaxed in the nuSUGRA 
scheme, where the non-universality in the Higgs sector allows for milder correlations among 
the parameters and among the low energy variables. This is especially true for the lower 
bound on tan (3, which now moves down to about 7, as displayed in Fig. 4b. Also for the 
other parameters we have weaker bounds with respect to the SUGRA case. For instance: 
niA < 450 GeV, mo > 200 GeV, and lightest stop and sbottom masses larger than about 
150 GeV and 400 GeV, respectively. 

In the effMSSM scheme most of the internal correlations of the model, which are typical 
of the supergravity-inspired schemes, are not present. For instance, tan/3 and itla are now 
independent parameters. Also in this case, the most relevant information about configura- 
tions which give Px and VLmh? inside the favoured region defined above is provided by the 
tan jS-rrih correlation, which is shown in Fig. 4c. We notice that in effMSSM the lower 
bound on tan/3 is around 5. In Fig. 4c we also display, by a dashed line, what would be 
the boundary of the scatter plot, in case set 2 for the quantities nig < qq >'s is employed 
instead of set 1. As for other relevant correlations, we find, for instance: < 500 GeV 
and/or rriq < 300 GeV. 

As a final comment, we point out that in our calculations we have taken into account the 
experimental constraint on sin^(a; — /?) This limit is stronger than the one displayed in 
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terms of tan (3 versus rrih in Figs. 4a-c and has the effect of depopulating the scatter plots, 
without modifying their boundaries. 

We wish to point out that, should the continuation of the LEP running provide some 
support in favor of a Higgs boson at a mass of about 115 GeV this would fit remarkably 
well within the configurations displayed in Figs. 4 [^]. It would be a case where accelerator 
measurements and WIMP searches would complement each other in providing a possible, 
consistent picture of physics beyond the standard model. 

For sake of comparison with specific experimental results, we provide in Figs. 5a-c the 
scatter plots for the quantity ^cricailr"^^ versus in the various supersymmetric schemes. 
^ is taken to be ^ = mm{l, Q^h"^ /{Q^h"^) mm}, in order to have rescaling in the neutralino 
local density, when flxh"^ turns out to be less than {flmh'^)min (here {flmh'^)min is set to 
the value 0.05). In universal SUGRA our results reach a maximum for ^CTgcaiar'"'^'' the 
level of about 10~^ nb, a feature which is in common, for instance, with the results of Refs. 
T3| , |29|j36| . Lower values for the WIMP-nucleon cross section are found in evaluations where 



various inputs, each one having the effect of suppressing the value of ^c'■ic^lar°'^^ ^^^^ employed 
concomitantly [|1^,0: i) low values for tan/5, tan/? < 10; ii) small values for the quantity 
rUs < ss >, iii) a tight lower bound on the neutralino relic abundance Q^h"^ > 0.1 . In fact, 
should we use the same inputs as in [|T^ , we would obtain the scatter plot which stays below 



the dashed line displayed in Fig. 5a, in agreement with the results of Ref. ||12|| . 



However, we point out that, in general, in universal SUGRA, evaluations by various 
authors differ in some features, for instance in the position of the maximum of ^Cgcaiar 
terms of m-^. This is likely to be due to the fact that this strict scheme is very sensitive to 
the specific ways in which various constraints (for instance, b — s + 7) are implemented in 
the calculations. 

In Figs. 5a-c the solid line denotes the frontier of the 3a annual-modulation region 
of Ref. 0, when only the uncertainties in pi and in the dispersion velocity of a Maxwell- 
Boltzmann distribution, but not the ones in other astrophysical quantities, are taken into 
account. As discussed in the Introduction, effects due to a possible bulk rotation of the 
dark halo or to an asymmetry in the WIMP velocity distribution would move this boundary 
towards higher values of m^. Our results in Figs. 5a-c show that the susy scatter plots reach 
up the annual-modulation region of Ref. [0], even with the current stringent bounds from 
accelerators (obviously, more easily in effMSSM than in SUGRA and nuSUGRA schemes). 

In connection with the results shown in Fig. 5c for effMSSM, we further remark that, 
if configurations with the hierarchy rriq < mi are included, the scatter plot would rise by a 
factor of a few at ~ 50 — 90 GeV. 

Finally, we recall that use of set 2 for the quantities < qq >'s instead of set 1 would 
entail an increase of about a factor 3 in all the scatter plots of Figs. 5. 



IV. CONCLUSIONS 

In this work we have shown that the current direct experiments for WIMPs, when in- 
terpreted in terms of relic neutralinos, are indeed probing regions of the supersymmetric 
parameter space compatible with all present bounds from accelerators. We have quantified 
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the extent of the exploration attainable by WIMP direct experiments in terms of vari- 
ous supersymmetric schemes, from a SUGRA scheme with unification assumptions at the 
grand unification scale to an effective model, effMSSM, at the electroweak scale. It has 
been stressed that, due the large uncertainties in the unification assumptions in SUGRA 
schemes, the effMSSM framework turns out to be the most convenient model for neutrahno 
phenomenology. 

Wc have proved that part of the configurations probed by WIMP experiments, and not 
disallowed by present accelerator bounds, entail relic neutralinos of cosmological interest. 
As discussed in the previous Section, this result is at variance with the conclusions of some 
analyses recently appeared in the literature. Also neutralinos which might contribute only 
partially to the required amount of dark matter in the universe have been included in our 
analysis. The cosmological properties have been displayed in terms of a plot of the local 
density versus the average relic abundance, i.e. in a representation which proves particularly 
useful to summarize the properties of relic neutralinos. 

We have noticed that a Higgs with a mass of about 115 GeV, such as the one now under 
experimental investigation at LEP2, would fit remarkably well in the above scenario. 

In our evaluations we have taken into account that the determination of the actual 
sensitivity region in terms of the WIMP-nucleon cross section and of the WIMP mass from 
the experimental data depends quite sizeably on uncertainties of various origins, mainly: i) 
possible effects due to a halo bulk rotation and/or to asymmetries in the WIMP velocities 
distribution, ii) significant uncertainties in the determination of Higgs-quark-quark and 
neutralino-quark-squark couplings. We have stressed that all these effects have to be taken 
properly into account, when conclusions about comparison of theory with experiments are 
drawn. 
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FIG. la. Scatter plot of cl^^iJ versus ^l^h"^ for universal SUGRA. Set 1 for the 
quantities < qq >'s is employed. Only configurations with positive /i are shown and 

is taken in the range of Eq. (|^). The two horizontal lines bracket the sensitivity region 
defined by Eq. (^. The two vertical lines denote the range 0.05 < Qmh"^ < 0.3. The region 
above Q^h"^ = 0.7 is excluded by current limits on the age of the universe. All points of this 
scatter plot denote gaugino configurations. 
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FIG. lb. Scatter plot of cr^cailr"^^ versus for nuSUGRA. Notations as in Fig. la, 

except that here the scatter plot contains neutrahnos of various configurations: dots denote 
gauginos, circles denote higgsinos and crosses denote mixed configurations. 
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Ic. Scatter plot of versus f^^^^ effMSSM. Notations as in Fig. 

;ns of fi are shown. 
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FIG. 2. Same as in Fig. Ic, except that here set 2 for the quantities < qq >'s is 
employed instead of set 1. 
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FIG. 3. Scatter plot of versus Q^h'^. This plot is derived from the experimental 

value [p^/{0.3 GeV cm~^) ■o'^^ailr°^^]expt = 1 ■ 10~^ nbarn and by taking in the range of 
Eq. (H), according to the procedure outlined in the text, in case of effMSSM. Set 1 for the 
quantities rriq < qq >'s is employed. The two horizontal lines delimit the range 0.2 GeV 
cm~^ < Px — ^-"^ GeV cm~^; the two vertical ones delimit the range 0.05 < ^Imh"^ ^ 0.3. 
The region above f^^^^ ~ ^-"^ excluded by current limits on the age of the universe. The 
band delimited by the two slanted dot-dashed lines and simply hatched is the region where 
rescaling of pi applies. Dots denote gauginos, circles denote higgsinos and crosses denote 
mixed configurations. 
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FIG. 4a. Scatter plot in the plane rrih — tan /3 of the SUGRA supersymmetric config- 
urations which stay either inside the region 0.2 GeV cm~^ < Px — GeV cm~^ and 
0.05 < Qmh'^ < 0.3 or within the corridor of rescaling in the plane versus n^h"^. Set 
1 for the quantities nig < qq >'s is employed. Crosses (dots) denote configurations with 



> 0.05 {n^h^ < 0.05). 



The hatched region on the right is excluded by theory. 

391. 



The 
The 



hatched region on the left is excluded by present data from LEP |^6|] and CDF 
solid line represents the 95% C.L. bound reachable at LEP2, in case of non discovery of a 
neutral Higgs boson. 
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FIG. 5a. Scatter plot of Co"scriar°°^ versus in case of universal SUGRA. Set 1 for the 
quantities rrig < gg >'s is employed. Crosses (dots) denote configurations with f^^/i^ > 0.05 
{Q^h'^ < 0.05). The dashed line delimits the upper frontier of the scatter plot, when the 
inputs of Ref. [jl2| are used. The solid contour denotes the 3cr annual-modulation region of 
Ref. Pi (with the specifications given in the text). 
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FIG. 5c. Same as in Fig. 5a in case of efFMSSM. 
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